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Abstract 



Consider a piecewise smooth expanding map of the interval possessing several in- 
variant subintervals and the same number of ergodic absolutely continuous invariant 
probability measures (ACIMs). After this system is perturbed to make the subinter- 
vals lose their invariance in such a way that there is a unique ACIM, we show how 
to approximate the diffusion coefficient for an observable of bounded variation by the 
diffusion coefficient of a related continuous time Markov chain. 
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1 Introduction 

Metastable dynamics arise in a number of physical systems. In such systems, the phase space 
can be divided into a finite number of components, called metastable states, that are nearly 
invariant under the dynamics. A typical trajectory will remain in one metastable state for 
an extended period of time before escaping to another metastable state and repeating this 
behavior. Rigourous results about metastability in dynamical systems perturbed by noise 
can be found in [B] and [9|. 
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In this paper we are concerned with reducing the description of purely deterministic 
chaotic systems to cooresponding finite state Markov chains. In particular, we continue 
the study of the dynamics of hyperbolic interval maps with metastable states initiated in 
[I]. These systems arise from perturbing an initial system To with m disjoint invariant 
intervals I±, I 2 , . . . , I m . The initial map has m mutually singular ergodic absolutely continuous 
invariant measures (ACIMs), fii,^, ■ ■ ■ , AW T is perturbed in such a way that the fij lose 
their invariance, and the perturbed map T e has only one ACIM, [i e . See Figure [TJ 



Such metastable systems can be understood in the context of deterministic dynamical 
systems with holes (see |3j) as follows. As the invariance of initially invariant intervals is 
destroyed by the perturbation, we think of the small set of points IiV\T~ x Ij that switch from 
Ij to Jj after the application of T e , as being holes in the initially invariant sets. Therefore 
the techniques developed to study systems with holes are useful in our analysis. 

As was shown in [7], we are able to approximate /i e , for small e, by a convex combination 
Y^j=iPj^3 °f the initially invariant measures, where (pi,P2, • • • ,p m ) is the invariant measure 
for the continuous time Markov chain on m states with transition rates proportional to the 
asymptotic sizes of the holes. Intuitively, since our map is chaotic on each interval, we 
expect the transition times between the Ij's to be almost independent, which explains the 
appearence of the above mentioned Markov chain. In this paper, we show that the Markov 
approximation also extends to the diffusion matrix for smooth observables. We believe the 
methods developed in our paper can be used to describe the transport coefficients in other 
chaotic systems, for example, billiards with narrow tunnels [12]. However, in order to present 
the ideas of our proof in the simplest possible setting, we restrict our attention here to the 
setup of [7J. 




Figure 1: A map with three almost invariant intervals 
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2 Statement of the main results 



In this section, we define a family of dynamical systems with m nearly invariant (metastable) 
subsets. They are perturbations of a one-dimensional piecewise smooth expanding map with 
m invariant sub intervals I\, I2, ■ ■ ■ , I m of positive Lebesgue measure. On each of these inter- 
vals, the unperturbed system has a unique ACIM. The perturbations break this invariance 
in such a way that each perturbed system will have only one ACIM. Our main result is an 
asymptotic formula for the diffusion coefficent of a smooth observable as the size of the per- 
turbation tends to zero. We also show that the sequence of jump times in between different 
intervals asymptotically approach exponential random variables. 

Let / = [0, 1]. In this paper, a map T : I is called a piecewise C 2 map with C = {0 = 
Co < c\ < ■ ■ ■ < Cd = 1} as a critical set if for each i, T|( CiiCi+1 ) extends to a C 2 function 
on a neighborhood of [q,q+i]. We call T uniformly expanding if its minimum expansion, 
inf;,; 6 /\eo I Tp (a;) I, is greater than 1. As is customary for piecewise smooth maps, we consider 
T to be bi-valued at points q G C where it is discontinuous. In such cases we let T(cj) be 
both values obtained as x approaches q from either side. 

2.1 The initial system and its perturbations 

The unperturbed system is a piecewise C 2 uniformly expanding map T : I with Co = {0 = 
c o,o < c i,o < ■ • ■ < Cdfi = 1} as a critical set. There are boundary points B = {bj} C (0, 1) 
such that Ij = [bj-\,bj] (b = 0, b m = 1) i.e. T (Jj) C Ij. The existence of an ACIM of 
bounded variation for T \j j is guaranteed by pi)]. We assume in addition: 

(I) Unique mixing ACIMs on the initially invariant subsets: To|/., j 6 {1 . . .m}, has only 
one ACIM /ij, whose density is denoted by <pj = dfij/dx. (T ,/ij) is mixing. 

From (I), it follows that all ACIMs of T are convex combinations of the ergodic ones, 

W- 

We define the points in Hq = (Tq 1 B) \ B to be infinitesimal holes. (The exclusion of 
boundary points from the set of infinitesimal holes is not essential, although it does simplify 
our presentation. See assumption (V) and the discussion thereafter.) 

(II) No return of the critical set to the infinitesimal holes: For every k > 0, (TqCq)P\Hq = 0. 
Since <pj are of bounded variation, they can be suitably defined so that they are continuous 

except on at most a countable set of points where it has jump discontinuities. Moreover (see 
[3 Section 4.2]), (II) implies that after <f>j have been so defined, they are continuous at each 
of the infinitesimal holes. 

(III) Positive ACIMs at infinitesimal holes: <pj is positive at each of the points in H fl Ij. 
(II) and (III) are generic conditions, although (II) may be difficult to verify for specific 

examples. Similar assumptions are made elsewhere, including in [91 Section 3.2] and [7j. 

(IV) Restriction on periodic critical points: Either 

(a) inf xeACo \T^(x)\ > 2, or 

(b) T has no periodic critical points, except possibly that or 1 may be fixed points. 

Because T may be bi-valued at points in Co, a critical point c^o is considered periodic if 
there exists n > such that c^o G T^jc^o}. Condition (IV) is necessary in order to ensure 
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that the transfer operators cooresponding to the perturbed systems defined below satisfy 
uniform Lasota-Yorke inequalities. These uniform inequalities are essential for establishing 
the perturbative spectral results, Propositions [2] and [3j which are a key ingredient of our 
proof. Since we cannot exclude the possibility of the forward orbit of a critical point contain- 
ing other critical points, these uniform inequalities do not follow directly from the original 
paper [TU], but rather from later works [21 H], see [7J Section 4.2]. 

In what follows, we consider C 2 -small perturbations T £ : / O of T for e > 0. A critical 
set for T £ may be chosen as C £ = {0 = c 0j£ < ci )£ < • • • < Cd, £ = 1}) where for each i, e \-t Ci j£ 
is a C 2 function for e > 0. Furthermore, there exists S > such that for all sufficiently small 
e, there exists a C 2 extension T i>£ : [c ifi - 5, c i+1>0 + 5] -> R of T £ \ [ci >e)Ci+li£ ] , and T ij£ ->■ T i)0 in 
the C 2 topology. 

We also assume: 

(V) The boundary points do not move, and no holes are created near them: Precisely, for 
each b G B we have the following 

(a) If b ^ Co, then necessarily T (b) = b. We assume further that for all e > 0, T e (b) = b. 

(b) If b G C , we assume that T (6_) < b < T (b + ), and also that b G C e for all e. 

This boundary condition can be considerably relaxed by suitably redefining the holes dis- 
cussed below, as was explained in [7J Section 2.4]. For simplicity of presentation, we do not 
detail these generalizations here. 

Set Hij j£ — Jj fl T e _1 (/j). We refer to these sets as holes. Once a T £ -orbit enters a hole, it 
leaves one of the invariant sets for To and continues in another. As e — > 0, the holes converge 
(in the Hausdorff metric) to the infinitesimal holes from which they arise. Our assumptions 
imply that there exist numbers > such that 

H i (H ij j S ) = ePij + o(e). (1) 

Conisider a continuous time Markov chain with states 1,2, ... ,m and jump rates from state 
i to state j equal to /?^. 

(VI) Irreducibility : The Markov chain defined above is irreducible. 

Condition (VI) implies that for small e > 0, T £ has only one ACIM fi £ , with density 
4> £ = d\x £ jdx. 

Examples of families T £ that satisfy (I) through (VI) can be found in [TJ Section 2.4]. 

It is natural to inquire if the family of ACIMs /i e has a unique limit as e — > 0, and if 
this limit exists, how to express it as a convex combination of /ii, . . . , // m . This problem was 
successfuly addressed in [7]. 

Proposition 1 (Theorem 1 in [7]). As e — )■ 0, 

(j) £ 00 = ^2pj<t>j 

j 

where (pi,p 2 , . . . ,p m ) is the invariant measure for the Markov chain. 
We define /xq to be the measure whose density is 4>q. 
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2.2 The main results 

Proposition [T] shows that the limiting Markov chain provides useful information about our 
metastable system. In this section we prove two additional results making the connection 
between the dynamical system and the Markov chain more precise. 

To state our first result we need to define, both for the T e -dynamics and for the Markov 
chain, a sequence of times that indicate when a transition occurs between different iy's. 
For the Markov chain: Set t^f = 0, and for % > 0, let tf* be the i th time that the Markov 
dynamics have changed states. Then for i > 1, set 7~ M = tf 1 —tfi v Let zf 4 be the state of the 
chain after the i th transition. Let P r denote the probability measure constructed on the space 
{1, . . . , m}' ' 00 ) for the Markov chain by starting at I r at t = and then evolving forward 
by the Markov dynamics. Then, given the present state of the chain, 7^ M are exponential 
random variables. That is, for t > 0, the densities are given by 

dP r (X M = t\z£ 1 =j) = p j e-e' t dt 

where (3j = J2k Pjk- Also given that = j, Z{ is independent of 1~ M and 

For T £ : Let z(x) = k if x e h- Set t§ = 0, and for % > 0, let tf = inf {n > t\_ x : z(T™x) ^ 

z{Te Ll x)}. Then for % > 1, set % £ = t\ - t\_ v 

Our first main result states that the finite dimensional distributions of jumps of the 
deterministic systems converge to the finite dimensional distributions of jumps of the Markov 
chain. 

Theorem 1. Fix j,p and S. For any intervals = [a*-, and numbers rjt £ {1, . . . , m}, 
k = 1,. . . ,p 

fijieTk 6 e A fc , z{t%) = r k , for k = 1, . . . ,p) -> F r {T k G A fc , z? = r k , for k = 1, . . . ,p) 
and the convergence is uniform for max b k < S. 

Consider an observable A : I — > R. For each fixed small e > 0, the Ergodic Theorem 
provides a law of large numbers, i.e. iV -1 ^2 k=0 A o T e fc — > fx £ (A) a.e. and in L 1 as N — > 
oo. If A is of bounded variation, the Central Limit Theorem applies [11] and states that 
N~ l l 2 *Y^ =Q A o T e fc approaches a normal distribution as N — > oo. We let D e (A) be the 
diffusion coefficent, which is the variance of the limiting normal distribution. 

Theorem 2. For any observable A of bounded variation, 

eD £ (A) -> D A/ (A) 

as e — > 0, where A is the observable on the state space of our Markov chain such that 
A(j) = Jj Acpjdx and D M stands for the diffusion coefficient. 
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Remark. D M can be computed efficiently. Assume for simplicity that A has zero mean 
(otherwise we subtract a constant from A). Let G denote the generator matrix of our 
Markov process, that is Gjk = (3jk if J '7^ k and Gjj = —f3j. Then 

D M = YV-A(j) / Pjk (t)A(k)dt = (pA, / e tG Adt) = ( P A, G^A). 

The proof of Theorems [T] and [2] depend on two perturbation results, Propositions [2] and 
[3j Proposition [3] is a variation on a result of Keller and Liverani [HI E] , see Appendix [A] 
for details, while Proposition [2] appears to be new. Our paper is organized as follows. In 
Section[3]we recall the transfer operator approach to the study of piecewise expanding interval 
maps and state Propositions [2] and |3j In Section [4] we derive Theorem [T] from Proposition [3j 
Section [5] contains the derivation of Theorem [2] from Theorem [T] and Proposition |2j Finally 
in Section [6] we explain how Proposition [2] follows from Theorem [TJ 



3 Preparatory material 
3.1 Function spaces and norms 

We use Leb to denote normalized Lebesgue measure on / and L 1 to denote the space of 
Lebesgue integrable functions on I. For / : I — > R, let \f\ L i = Jj\f(x)\ dx, \f\ Loo = 
sup x&1 \f(x)\, and var(/) be the total variation of / over /; that is, 

n 

var(/) = sup{^ \f{xi) - f{xi-x)\ : n > 1, < x < xi < ■ ■ ■ < x n < 1}. 

i=i 

BV = BV(J) is the Banach space of functions /:/—>■ R with norm 

11/11 = t i nf „ var(^) + \g\ Lao . 

g=j except on a countable set 

An element of BV is technically an equivalence class of functions, any two of which agree 
except on a countable set; we generally ignore this distinction. 



3.2 Transfer operators and the spectral setting 

For e > 0, let C e be the transfer operator associated with T £ acting on BV, i.e. 

My) 



C £ A(x)= 



T'M\ 



Note that £0 has 1 as an isolated eigenvalue of multiplicity m, and then a spectral gap. Let 
P denote the associate spectral projection, that is 



(PA)(x) 



J A(y) d v \ ^jix) if x e h 



For small e > 0, 1 becomes an isolated eigenvalue creating a small spectral gap. Note that 
C £ preserves the space BVq of BV-functions with zero mean. 
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Proposition 2. There exists 77 < 1 and k > such that 



||^ e ||BVo S ^ • 

The proof of Proposition [2] is given is Section [6] 

Remark. In case of two intervals much more precise asymptotics of the spectral gap can be 
deduced from the results of [5]. It is likely that a similar statement holds for an arbitrary 
number of intervals, however we do not pursue this question here since the weaker version 
stated above is sufficient for our purposes. In fact we hope that the argument used to 
prove Proposition [2] can be extended also to the case of piecewise hyperbolic systems such 
as billiards consisting of several regions connected by small holes. 

For j G {1, . . . , m}, let Cj >£ be the transfer operator acting on BV(Ij) where we consider 
Hj £ as a hole, i.e. 

yeiijnTr 1 ^}) 



m(y)\' 



From [5] we see that, for small e > 0, Cj >£ has an isolated simple eigenvalue Xj >£ < 1 of 
multiplicity one with A J)£ — > 1 as e — > 0, and otherwise the spectrum has a uniform spectral 
gap. More precisely, if Qj >£ is the spectral projection corresponding to A Jj£ , the following 
statement holds. 

Proposition 3. There exists a measure Uj j£ and a function J)£ G BV(ij) such that 

VjMhe) = 1 , Qj,s = VjAWte 

and 

(a) X j)S = 1 - 0j£ + o(e). 

(b) As £ — > 0, Uj )B (A) — > f A(y)dy where the convergence holds in strong topology in BV* 
and (j)j j£ — > 4>j in L 1 . Moreover 

limsup 1 4>j,e - 07 1 = 0. (2) 

(c) There exists C > such that, for all small e, A G BV(Jj), and all n > 0, 

||A^ £ A-Q i)e A||<Cr||A||. 

This proposition follows from the work of Keller and Liverani [HI E] , see Appendix [A] for 
details. 
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4 Convergence of the jump process 
4.1 Tightness 

We need to know that the distribution of jumps is tight so we begin with the following result. 

Lemma 1. Given S, 5,p there exists a such that, for all e sufncently small, 

/i p (3Jfe with t% < S/e and T k £ +1 < a/e) < 5. 

The following estimate plays a key role in our analysis. Given a segment J G I let r n (x) 
be the distance of T™x to the boundary of the component of T™ J containing it. Recall that 
Leb denotes Lebesgue measure on / = [0,1]. 

Lemma 2 (Growth Lemma), (see [3j Section 5.10]) There exists A > 1, c > such that, for 
all e small enough, and all J, and all n > 0, 

Leb(x : r n (x) < e) < Leb(x : r (x) < e/A n ) + cLeb(J)e. 

Proof of Lemma^ First, t$ = 0, and if N = N(a, e) = \cr/s \ , then 

/i P (7^ < k/s) < 1 - /i p (Ti > N) — 1 — J dx = l-\^ £ J Q L , £ (0 L ) dx + O(0 N ), 

where we have used Proposition [3] (c). But from parts (a) and (b) of the same proposition, 
we see that this can be made arbitrarily small for all small e by taking o small enough. 

Because \i\ Leb it suffices to prove the statement for Leb. 

We follow [5], Section 18. Let S n<m = Y^=n+i ^-t 3 x^h e - We have to show that 




( x ) >0 dx = o(l), e ->■ 0, o ->■ 0. 



Take a small r. We say that a visit of x to the hole at time n is (r-)inessential if the length 
of the smoothness component of T"L fl H £ containing T £ x is less than re. By Lemma [2] the 
probability that x will have an inessential visit to the hole before time S/e is less than CrS 
which can be made as small as we wish by taking r small. Therefore it suffices to show that 
for any fixed r 

/ 1 ^,eW 1 S„, ff/£ (l)>0^ = o(l) 
n<S/e J 

where S n , £ (x) = {x has essential visit to the hole at time n}. 

J 1 ft,, s (a>) 1 5 Bjff/ ,(*)>oda; = Leb(£ n>£ (x))F(S nta/£ (x) > 0\£ n>£ (x)). 

Since Leb(S n ^ e (x)) < Leb(x : T™x G H £ ) < Ce it suffices to check that 

max Leb(S ntCT / £ > Q\£ nt£ (x)) = o(l), e — > 0, a — > 0. (3) 

n<S/e 
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In order to prove this we observe that due to assumption (II) for any fixed M , S nj M ( x ) = 
for all x G £ n .e provided that e is small enough. On the other hand if k > Mq then by 
Lemma [2] applied to H e 

Leb(l H T? +k x) > 0\(S n;£ (x)) <C(e + A~ k ) . 
Summing over k G [M + 1, k/e] we obtain □ 



4.2 Convergence of the finite dimensional distributions 

Proof of Theorem^ The proof is by induction on p. First, for p = 1 assume that the initial 
distribution of x is chosen according to some density p G BV(pj). Then by Proposition |3]^c) 

p J (T 1 £ = n,z(T?x)=r)= f C^(p)dx = X% f (Q 3 , £ p)dx + 0(6 n ). 
Recall that 

By Proposition |3j z/ Ji£ (p) — > J p(x)dx, J H 4>j, £ /e — > /3j r and for en w t, A™,, —> e~^ t . Thus 
summation over n G A/e concludes the proof for p = 1. 

Next suppose that the statement is known for some p. Denote 

Q = {T k £ G A k ,z{t £ k ) = r k , for k = 1, . .. 

To carry the induction step it is enough to prove that 



/ii(e7^i G A p+1 2;(£; +1 ) = r p+1 , fi) -> fij(n)^ ZpZp+1 / e 

•/A p+ i 



(4) 



Set 



(C n A)(x)= 



My) 



T £ e y=x 



where the sum is taken over y G fl. Then 

h((^V +1 = nz(t; +l ) = r p+1 , n) = f q^cuh^dx. 

Take a small a and rewrite 



(5) 



lz p z p+1 



Due to Lasota-Yorke inequality 
has bounded BV-norm. Hence arguing as in the p — 1 case we see that ^ is asymptotic to 



-p*ptdt 



•p+1 



£;/ £ (£n(l/,^)) rfx(l + o^ (l))- 



Due to Lemma [T] the last integral here equals to 

Since <r is arbitrary this proves Q completing the proof of Theorem [I] 



□ 
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5 Proof of Theorem 



In this section we deduce Theorem [2] from Theorem [T] and Proposition |2j 
Proof. Let A = A — ft £ (A). By Proposition for each 5 we can find S such that 



E MAAoT £ 

n>S/e 



< 5 



so it suffices to get the asymptotics of /j, £ (AA o T") for n ~ t/e where t < S. We need to 
estimate D n = fi £ (AA o T"). Take a natural number no- We have 



D n = j A(x)4> £ (x)A(T™x)dx = j A(T?y) [C^CCW.)] (v)<*V- 
Since L £ depends continuously on e as a map BV — > L 1 we can rewrite the last expression as 

A(T?°y) [q:- na {P{A4> £ ))] (y)dy + O(9 no ) + o e ^ (l). 
Since To is mixing we have 

J A(T?°y)B(y)dy = J2( K J i Mv)<f>kdy jf + [O(0"°) + o e _*(l)] |L4| | BV | \B\ | BV - 

As e — >• the first factor converges to A(fc) — ^2iPiA-{l) while the second term equals 



n - 2n ^ dy + o £ ^{\) 



Since n was arbitrary we conclude that 



w e ^(^°^)=e r 



Pj A(j)A(k) Pjk (t) - Eft A 0') 



n<5/e jfc 

Since S 1 is arbitrary we can let 5 — > oo and obtain Theorem [2] 



□ 



6 Proof of Proposition [2 

Proof. Recall that £ £ satisfy uniform Lasota-Yorke Inequality, that is there is a constant K 
such that 



var(£™A) < K [\\A\\ L , + ^var(A)] 
||CA|| L i <K\\A\\ L , 



(6) 
(7) 
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Consequently it suffices to show that if A e BV then 

\\£ K J 2e A\\ Ll < (^^PHbv (8) 
since combining (J6]), ([7]) and rt8J) gives 

\\£ K £ / £ A\\ m <2- 1 \\A\\ m . 
Next mixing of To on implies that if Jj A(x)dx = for each k then 

\m°A\\ L i<ce no \\A\\ BV . 

Consequently given 5 we can find n ,e and 5\ such that if e < e and | f T A(x)dx\ < 5i for 
each k then 

\\£?\\v<5. 



Let n = Y e ~~ n o- By the foregoing discussion it remains to show that for each k 



[C n e A)(x)dx 



<8x 



provided that k is large enough. Since 

C n £ A = C n - ni C r ^A = C^ ni {PA) + 0{6 ni ) + o e ^ (l) 

we need to show that J /fe (£™~ ni PA)(x)dx is small. By Theorem jlj this integral is asymptotic 
to ^2 ■ fj A(y)dypjk(K). As k — > oo this expression converges to L A{y)dy = so the result 
follows. □ 



A Proof of Proposition [3 

Proof. Part (a) is proven using arguments similar to the ones in [9} Section 3.2]. Namely we 
apply P, Theorem 2.1] with Vo = Co, V e = £j, e , and the Banach space BV(Jj). This theorem 
says that 



lim - s 



1 - > ,Qk 

k>0 



with 



A e = idj(H j)e ) = fye + o(e), 

- lim ''^ 1 ) - ^ {k+l)R ^ 

Now for fixed k and all small e the set T~ k Hj )£ consists of a finite number of intervals of 
size 0(e) so the sizes of their preimages by T and T e differ by 0(e 2 ) (here we are using 
assumption (II), which implies that 0o is continuous at points in T H £ ). 

This completes the proof of part (a). Parts (b) and (c) follow from [8] except for ^ 
which is proven in [7]. □ 
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